@sdat2 , thanks for your contribution !

this new function seems to return an unravelled version of the preprocessed data, which in this specific PCM configuration, happens to be the stacked/joined PCA eigen values.

so, my comments:

1. if users modified default options, the returned array may not corresponds to the PCA value (eg: if using another reduce function), so the vocab here would be wrong.
2. the pca dimension of the output array could be generalized as a "reduced_dimension", which is a reduced version of the dim passed as argument.
3. I would rather return by default a dataset where the feature variables are back into separate dataArrays, instead of the stacked version here. At least this should be an option (and possibly takes into account the tentative future PCM option "join" discussed in #27 )

The request changes are thus a more general version of this function taking into account the 3 comments above.
Then the documentation could be updated to show users how to get access to PCA values using this new function.